Global existence and future asymptotic behaviour 
for solutions of the Einstein- Vlasov-scalar field 
system with surface symmetry 



We prove in the cases of plane and hyperbolic symmetries a global 
in time existence result in the future for cosmological solutions of the 
Einstein- Vlasov-scalar field system, with the sources generated by a dis- 
tribution function and a scalar field, subject to the Vlasov and wave equa- 
tions respectively. The spacetime is future geodesically complete in the 
special case of plane symmetry with only a scalar field. Causal geodesies 
are also shown to be future complete for general homogeneous solutions 
of the Einstein- Vlasov-scalar field system with plane and hyperbolic sym- 
metry. 

1 Introduction 

Consider the Einstein- Vlasov-scalar field system with spherical, plane and hy- 
perbolic symmetries. For more information on this system which describes the 
evolution of self-gravitating coUisionless matter and scalar waves within the con- 
text of general relativity, see We study in this paper global existence for 
solutions of this system where the Einstein equations are coupled to the Vlasov 
and wave equations. This completes in the expanding direction, the results of [H] 
which concentrated on the local in time existence results and also in the global 
existence in the special case of plane symmetric Einstein-scalar field system (i.e 
without Vlasov contribution). It is shown that in the spherical case, no solution 
exist globally in the future. Using the results of 3 which studied the future 
asymptotics of polarized Gowdy T"^-models, we prove in the plane symmetric 
case without Vlasov that the spacetime obtained is geodesically complete in the 
future. The same result holds for spatially homogeneous solutions of the full 
system with plane or hyperbolic symmetry. Only a linear scalar field is treated 
but it is likely that the approach developed here will be useful in the nonlinear 
case. 
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Abstract 



Let us recall the formulation of the Einstein- Vlasov-scalar field system as 
shown in 9 . We consider a four-dimensional spacetime manifold M, with local 
coordinates {x") = {t, x*) on which x'^ = t denotes the time and (a;*) the space 
coordinates. Greek indices always run from to 3, and Latin ones from 1 
to 3. On M, a Lorentzian metric g is given with signature (— , -f , We 
consider a self-gravitating coUisionless gas and restrict ourselves to the case 
where all particles have the same rest mass, normalized to 1, and move forward 
in time. We denote by (p") the momenta of the particles. The conservation 
of the quantity gapp'^P^ requires that the phase space of the particle is the 
seven-dimensional submanifold 

PM = {5a/3pV = -1; / > 0} 

of TM which is coordinatized by (i, If the coordinates are such that 

the components goi vanish then the component is expressed by the other 
coordinates via 

p" = ^^^Ji + gvP'p' 

The distribution function of the particles is a non-negative real- valued function 
denoted by /, that is defined on PM. In addition we consider a massless scalar 
field </> which is a real- valued function on AI. The Einstein- Vlasov-scalar field 
system now reads: 

V"Vc,(/) = 

Gap — STrTap 

f 1 dp^dp^dp^ 1 

Tap = ~ fPaPp I 5 I " ^ (Va(/>V,3(/) - -ffa^gV^C^Vc/)) 

where Pa = gapP^ , \g\ denotes the modulus of determinant of the metric gap, 
T^p the Christoffel symbols, Gap the Einstein tensor, and Tap the energy- 
momentum tensor. 

Note that since the contribution of / to the energy-momentum tensor is 
divergence-free [2], the form of the contribution of the scalar field to the energy- 
momentum tensor determines the field equation for (j). 

We refer to |7j for the notion of spherical, plane and hyperbolic symmetry. 
We now consider a solution of the Einstein- Vlasov-scalar field system where all 
unknowns are invariant under one of these symmetries. We write the system in 
areal coordinates, i.e coordinates are chosen such that R = t, where R is the 
area radius function on a surface of symmetry. The circumstances under which 
coordinates of this type exist are discussed in HI . In such coordinates the metric 
g takes the form 

ds^ = + e^^f^'^-^dr^ + t^dO^ + sin^. Od^^) (1.1) 
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where 

{sin 6 for k = 1 (spherical symmetry); 

1 for fc = (plane symmetry); 

smh0 for fc = — 1 (hyperbolic symmetry) 

t > denotes a time-like coordinate, r e M and {9, ip) range in the domains 
[0, tt] X [0, 27r], [0, 2tt] x [0, 2tt], [0, oo[x [0, 27r] respectively, and stand for angular 
coordinates. The functions A and /i are periodic in r with period 1. It has been 
shown in |S] that due to the symmetry, / can be written as a function of 

t, r, w := e^pi and F := t^lip^f + sin^ 9(j)^f], 



i.e. / — f{t,r,w,F). In these variables, we have p° = e^^\/l + w"^ + F/t^. 
The scalar field is a function of t and r which is periodic in r with period 1. 

We denote by a dot and by a prime the derivatives of the metric components 
and of the scalar field with respect to t and r respectively. Using the results of 
P], the complete Einstein- Vlasov-scalar field system can be written as follows: 



a/ 1 -f- ty^ + /i^ 

e-^''{2tX + l) + k = 8nt^p (1.3) 

e"2^(2t/i- 1) - fc = STTt^p (1.4) 

fi' = -47rie^+^j (1.5) 

e-2^M" + Ai'(M' - A')) ~ e-2^(A + (A + i)(A - A)) - 4nq (1.6) 

e-2A^" „ e-2M^ „ _ ^ + ^)0 _ e'^^iX' - = (1.7) 
where (jl.Tfl is the wave equation in and : 



OO /' + 00 



p(i, r) = e-2'^roo(t, r) = ^ I / ^1 + + ^/^^/(^, ^, F)di^du; 



+ i(e-2M^2^e-2V'') 



(1.8) 



J-oo Jo ^1 + w2 + F/t^ 



(1.9) 
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(1.10) 



TT f°° f°° F 



q{t,r) = j^T22{t,r) = -^—^T33it,r,e) 



, , , , =f(t,r,w,F)dFdw + e-^''<p^ ~ e-^^(f>'^ 

(1.11) 

We prescribe initial data at time t — \: 

fil,r,w,F)^f{r,w,F), A(l,r)=A(r), /i(l, r) = ^(r), 

O 

0(1, r) = (t>{r), r) = i;{r) 

The choice < = 1 is raade only for convenience. Analogous results hold in the 
case of prescribed data on any hypersurface t ~ to > 0. 

The paper is organized as follows. In section 2, we show that the solution 
of the Cauchy problem corresponding to system H1.2() - (|l.ll|l exists for alH > 1 
and fc < 0. In section 3, we prove that the spacetime obtained is timelike and 
null geodesically complete towards the future in the plane symmetric case with / 
identically zero. This section ends with results on the homogeneous models. We 
prove in section 4, for plane symmetric solutions with / = 0, that the maximal 
globally hyperbolic development of any initial data is inextendible. 

Unless otherwise specified in what follows, constants denoted by C will be 
positive, may depend on the initial data and may change their value from line 
to line. 

2 Global existence in the future 

In this section we make use of the continuation criterion in the following local 
existence result. 

Theorem 2.1 Let / G C\R'^ x [0,oo[) with 

}{r + 1, w, F) = /(r, w, F) for (r, w, F) G x [0, c»[, / > 0, and 

O 

Wo :— snp{\w\\{r, w, F) G supp/} < cx) 

o 

Fo := sup{F|(r, w, F) e supp/} < oo 
Let A,V' e C^R), C2(R) with A(r) = A(r + 1), /i(r) = ^^(r + 1), 

o o 

(j){r) = (j){r + 1) and 

ll{r) = -47re^+^}(r), r G R 
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Then there exists a unique, right maximal, regular solution (/, A, fi, (p) of system 



wUh (/,A,Ai,0)(l) = {f,X,fi,(j)) and 0(1) = V 

on a time interval 

[1,T[ with T e [l,oo]. // 



sup{| w I \{t,r,w,F) e supp/} < oo 



and 



then T — oo. 



sup{e 



2fj.(t,r) 



e M, t e [i,r[} < cx^ 



This is the content of theorems 4.6 and 4.7 in 0. For a regular solution, all 
derivatives which appear in the system exist and are continuous by definition 
(see |5]). We prove the following result in the cases of plane and hyperbolic 
symmetries. 

Proposition 2.2 Assume that (/, A, /z, 4>) is a right maximal regular solution of 
the Einstein- Vlasov- scalar field system with initial data given for t — 1. Then 

for k — or k = —\, 



sup{e 



2fj.{t,r) 



|r G R, te [1,T[} < oo. 



Proof : We now establish a series of estimates which will result in an upper 
bound on /j. Similar estimates were used in ;9 and 1\. Firstly, integration of 
(|1.4|l with respect to t and the fact that p is non-negative imply that 



t 



t 



s^p{s, r)ds 



> 



t 



Co - kt 



te[i,T[ 



(2.1) 



where Co ^ k + sup{e~^^('''; r £ [0, 1]}. Next let us show that 

[ e''+^p{t,r)dr <Ct, t £ [1,T[ (2.2) 
Jq 

A lengthy calculation shows that, since fi and j are periodic with respect to r : 



2p + q-^{l + ke'n 



dr (2.3) 



For fc = 0, since q > — 2p and p < p, we have: 



4/ et'+^p{t,r)dr <- [ e''+^ p{t,r)dr 
at Jo t Jo 
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Integrating with respect to t yields (|2.2|l for fc = 0. For fc = — 1, we have, using 
1)2. l|l and the fact that 2p + q > — p, p < p: 



dt 



e^+^pit, r)dr 



2p + q-P±^{l-e^n 



dr 



< - 



2 ^ C'o + t Jq 



e^+^P±Pdr 
2 



<( 



t Co + t Jo 



) / e^'+^pdr 



Integrating with respect to t yields (|2.2|l for fc = — 1. Using 1)1. 5|) . the fact that 
|j| < p and 12.2|l we find 



I ^iit,r) - / p{t,CT)da I =1 



1 pi 



Ja 



T)\dTda 



Jo 



that is 



< Ant [ e>'+^p{t,T)dT 
Jo 

I p{t, r)- [ pit, a)da \< Ct^ , t e [1, T[, r G [0, 1] 
Jo 



Next we show that 

By H2.1|l and since p > p, 



(2.4) 



(2.5) 



at 



ATTte^f^ip- p) 



1 + fce^^ 



t 



<(t + 



1 



t Co-kt' 



and integrating this inequality with respect to t yields 



sf"^ < C 



t 



Cn - kt 



< Ct, 



i.e. (ESJ. 

Now we estimate the average of p{t) over the interval [0, 1] which in combi- 
nation with (|2.4(l will yield the desired upper bound on p. We use H1.4II . I|2.2|l . 
ifT^ and the fact that p < p, fce^'^ < : 



t /.I 



p(t,r)dr — i p{r)dr + 



/i(s, r)drds 



< / p{r)dr + - In i + 47r 
Jo 2 

/•I 1 

< / °p{r)dr + -\nt + Ct^ 
Jo 2 



se'^-^e^'+^pdrds 



1 Jo 



1 JO 



fce^'^ 



drds 



6 



with (|2.4|l this implies 

/i(i,r) < Cil + lTLt + f^ + t^) < Ct^, t e [1,T[, r e [0,1]. (2.6) 



Remark 2.3 1/Fe /lawe proven that for initial data as in theorem \2.1[ the right 
maximal regular solution of the full system satisfies estimates 12.^) . \2. 5]) and 

In the next theorem, we prove that the solution exists on the full interval [1, c»[. 

Theorem 2.4 Assume that {f,X,fi,(f)) is a solution of the full system on a 
right maximal interval of existence [1,T[ with initial data as in theorem \2.U 
then T = oo. 

Proof: Assume that T < oo. We show that under the bound of /i, we obtain the 
bound of suplluij r, w, _F) G supp/} which is in contradiction to theorem 
12.11 The proof of the bound on w is similar to the proof of [see 0, theorem 
6.2]. 

Except in the vacuum case we have wq > and Fq > 0. For t > 1 define 
P+{t) max{0, maxjuiKr, ui, F) e supp/(t)}}, 



P-{t) := min{0, minjuiKr, ui, i^) S supp/(t)}}. 

Let (r(s), w(s), F) be a characteristic in the support of /. Assume that P+{t) > 
for some t e [1,T[, and let w{t) > 0. We have 

w{s) = -Aw - e^'-^n'y/l + w^+F/s^ 
47r2 



„2m 



oo />CXD 



OO Jo 



2s 



-vu ■ 



1 



F/s^ -wJl + w^ + F/s^ fdFdw 



-w{(f) e 



-2/j 



h'^e-^^) + ^l + w-^+F/s'^^(j}'e-^'-^ 



(2.7) 
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Set 7 — w\/l + w'^ + F/ — wy 1 + w'^ + F / s^. As long as ti;(s) > 0, we have 
the foUowing estimates: if w < then 7 < 0. If > then 

47r2 



„2p 



00 nOO 



H 

s 

< e 

s 



„2m 



2^^ 



-00 JO 



^ fat aw H w < w 

n(-) rFo ^2(1 ^ ^,2 ^ p/p-^ -w^l + w^+ F/s^) 



fdFdw 



wy/T+w^^^F/s^ + wyl + w'^ + F/s 

P+{s) f^" w{l + F) ^ -^^ 1 + 

' jai'dw-i w 



2s 



o2m 



1 1 + kp'^^^ 
<in-^Fo{l + Fo)\\f\\—iP+is)f- + 



2s 



2s 



C 
< — 

s 



(2.8) 



It is proved in ^ that K{s) = sup{(e"2^02 +e"2^(^'2)(s, r); s G [1,T[} < K{l)s, 
then 



> (w + y^l + w^ + F/.i 



> -v/i + w2 + i^/s2iv:(s) 

> -C(l + P+(s))s 
Therefore, using ijTHl) and (jTH), ||22Il yields: 



(2.9) 



w{s) < 



C 



w{s) 



- w{s) 



+ C(l + P+(s))s 



I.e. 



-^w{sf < C(s-i + s)(P+(s))2 + CP+{s)s 
ds 

as long as w{s) > 0. Let ti G [1, i[ be defined minimal such that w{s) > for 
s G [ti, t[, then 

witf < wihf + C f\is-' + s)(P+(s))2 + sP+is)]ds 

If ti = 1 then w{ti) < wq and 

w{tf <wl + C [(s-i + s)(P+(s))2 + sP+{s)]ds. 
If ti > 1 then w{ti) = and 

witf <C f [(s-i + s)(P+(s))2 + sP+{s)]ds 



<wl + C / [{s-^ + s){P+{syf + sP+{s)]ds. 
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Thus 

(P+(t))2 <wl + C j\{s-' + s){P+{s)f + sP+{s)]ds 
<wl + C [(s-i + .s)(P+(.s))2 + s]ds 

< {wl + Ct^) + C {s-^ + s){P+{s)fds; for all t G [l,r[ 
Applying Gronwall's inequality to this estimate implies that P+ is bounded on 

Estimating w{s) from below in the case w{s) < along the same lines shows 
that P_ is bounded as well and the proof is complete. 

Remark 2.5 In the spherically symmetric case (k — 1), there is no global 
existence in the future, regardless of the size of initial data. For any solution 
(f,X,IJ,,(j)), the estimate 

e-2''(*''^) = ^— i - 1 - y y s'v{s,r)ds < ^—^ - 1 

has to hold on the interval [1,T[. Since the right hand side of this inequality 
tends to —1 for t —i- oo, it follows that T < oo. Furthermore we deduce from the 
previous theorem that ||e^^|| — + oo for t ~* T . 

3 The future asymptotic behaviour 

In this section we prove that the spacetime obtained in theorem l2.4l (with / = 0) 
in the plane symmetric case is timelike and null geodesically complete in the 
expanding direction. Later on, we prove that this result holds for homogeneous 
solutions of the Einstein- Vlasov-scalar field system in plane and hyperbolic sym- 
metries. 



3.1 Particular explicit solutions 

Let us determine first the explicit solutions <f>, fi, A in the case / = and k = 0. 
In the case k = the wave equation can be reduced to a simple linear 
equation and the result follows from [S]. This reduction goes as follows. In that 
case the field equations imply that X — ^ + logt is constant in time. It may, 
however, be dependent on r. Suppose that r is replaced by a new coordinate 
s on the initial hypersurface. Choosing s appropriately makes the transformed 
quantity X — ii + log t constant on the initial hypersurface and hence everywhere. 
Once this transformation has been carried out the wave equation simplifies to 0-1- 
t^^(f) — (f)" . We can use the results of Jurke [2j obtained in the case of polarized 
Gowdy T'^-models (where W is replaced here by 0, t by a multiple of t^ and it is 
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necessary to be careful to get the right periodicity of the transformed solution) . 
The equations are similar with the same boundary conditions. Following this, 
for all {t,r) G [l,oo[xR, (j) must be in this form : 



f «! + 26 logt (homogeneous case) 

I Oi + 26 logt + t^^a{t, r) + P{t, r) (non-homogeneous case) 

(3.1) 

where Oi and h are constants fixed by initial values of (p and 0, a and (3 real- 
valued C^-functions with |a|, \a'\ < C, |d| < Ct, \/3'\ < Ct'^, |/3| < 

and jit'^dtt - t~^dt)a{t,r) = drra{t,r). With this, < Ct-^. 
From the field equations (|1.3I) . 11.4|l we have : 

A = -^ + 27rt(02+4t20'2) and (i = + 2'Kt{^'^ + U"^ cf)'"^) 

Therefore /t and A are bounded each by a positive constant C . From Theorem 
15 of P] (the hypothesis of this theorem are satisfied : at is replaced here by /t), 
/X can be cast for all {t,r) e [1, oo[xR into the form 

I vt -f o(t, r) (non- homogeneous case) 

where v \s a, positive constant, 7 a constant fixed by initial value of /i and the 
function 5 satisfies the inequalities, |<5(t, J^)! < C(l +t), \5\ < Ct. Note 
that if (/()' = then from equation (1.5), ^' = and A' = 0; the solutions are 
independent of r and we are in the homogeneous case. 

3.2 Geodesic completeness 

Let ]t_, r+[9 T h-> (x"(t),p'^(t)) be a geodesic whose existence interval is max- 
imally extended and such that x°(ro) = t(To) = 1 for some tq e]T_,r+[. We 
want to show that for future-directed timelike and null geodesies, t+ ~ -\-oo. 
Consider first the case of a timelike geodesic, i.e., 

= -m^ ; / > 

with TO > 0. Since dt/dr = p° > 0, the geodesic can be parameterized by the 
coordinate time t. Recall that along the geodesies the variables 
t, r, p", w := e'^p^, F := [(p^)^ + sin^ ^(p'^)^] satisfy the following system of 
differential equations : 

dr X dw ; n 9u \ 1^ n\9 dF 



e W 



. ^, A/w-e2M-A^.(pO)2^ " ^0 (3 3) 

dr dr dr 

— =p°^ ^ = -/i(/)2 - 2e-Vy^ - e"'''Au;2 _ e-^'^t'^F (3.4) 
ctT dr 

With respect to coordinate time the geodesic exists on the interval [l,oo[ since 
on bounded ^-intervals the Christoffel symbols are bounded and the right hand 
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sides of the geodesic equations written in coordinate time are linearly bounded 
in p^, , p'^ . Along the geodesic we define w and F as above. Then the relation 
between coordinate time and proper time along the geodesic is given by 

dr 

and to control this we need to control w as a function of coordinate time. 



3.2.1 The plane symmetric case without Vlasov 



Assume that w{t) > for some t > 1. By the fact that |0||(/)'|(s) < Cs~\ 

g(M-^)(s) — (jg aj^c]^ |j| < have as long as 'w{s) > 

w{s) — 47rse^^(j-\/ rn^ + w'^ + F / — pw) + —w 



ds 



2s 



2s 



(3.5) 



< —w + Anse^'^^lSU'lJm^ + F/s^ 
2s 

1 

< —w + Cs 
2s 

Let to e be defined minimal such that w{s) > for s E [ta,t[. Then 

Gronwall's inequality shows that 



w{t) < 



w{to) + C sexp( 



to 



27 



)ds 



exp( 



dr , 
2^' 



Now either to — 1 and w{to) — 'w{l) or to > 1 and u'(to) = 0. Thus 



\wil)\+C I sexp(/ ^)ds 



Estimating w{s) from below in the case w{s) < along the same lines yields 
the upper bound of —w{t). 



Since < C{l + t), ut^ + S > vt^{l-^^^). For t large, we choose 



C(l+t) 



< 



Then along the geodesic we have 
dr 



> 



dt y/m^ + + F/t^ ~ 2tVm2 + C + i^ 

Since the left hand side is constant, the integral over [l,oo[ diverges. 
In the homogeneous case = 0, j — 0, p > 0), becomes 

d , ^ 1 
—w(s) < — w 
ds 2s 

i.e w{t) < foj. ^ > I Therefore 

dt ~ 
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The integral of the right hand side of the above inequahty over [1, oo[ diverges. 
In either case, we conclude that t+ = +00 as desired. 

In the case of a future-directed null geodesic, i.e. m = and p"(to) > 0, p° is 
everywhere positive and the quantity F is again conserved. The argument can 
now be carried out exactly as in the timelike case, implying that t_|_ = +00. We 
have proven : 

Theorem 3.1 Consider initial data with plane symmetry for the Einstein- scalar 
field system written in areal coordinates. Then the corresponding spacetime is 
timelike and null geodesically complete in the expanding direction. 



3.2.2 The spatially homogeneous case 

Consider now the Einstein- Vlasov-scalar field system for fc < 0. Wc study the 
future behaviour for solutions which are independent of the space coordinate r. 
In this case /i' = 0, </>' = and equation l|1.5|) shows that j — 0. 

We prove that w{t) is bounded for finite time t and the spacetime obtained 
is geodesically complete. Since p > 0, fc < 0, we have from 1)3. 3() . 

—w[s) — —Anse ^ pw H w < — w 

ds 2s 2s 

and we can prove as above that w{t) < Ct^^^ for t > 1. Using (|2.1|) . along the 
geodesic. 



dr 
It 



> 



C 



Ct + F/i2 



The integral of the right hand side over [l,oo[ diverges. Therefore r+ = -l-oo. 
In the case of a future-directed null geodesic, the argument is the same as what 
we have done in the previous subsection. 

In the homogeneous number of further estimates can also be obtained. 

We follow the proof of theorem l2.1l and 6 to obtain a better bound of \i which 
depend explicitly of the data. Recall (|2.3() 



2p 



p + p 



(1 + fce^f ) 



Since <; > 0, we obtain 
for fc = 0, 



thus 



e^+^p{t) < e^+^°pt-^- 



(3.6) 
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2p 



p + p 1 P + P 



2 Co + t 2 
1 1 



< 



Ci+t 



t Co + t 



where Co = k + e and Ci = max(0, Co). Thus 

e^+^p(t) < (Ci + l)2e^+^pt-2. 
Next we have similarly to (|2.5|l . 



e-2A' - 1 + t 



< el^l"^ for k = -1 



For fc = 0, 



p{t) < A* + ^ log i + 47r ^ se^^'pds ^ p + ^\ogt + Aw 



< ^ + - log i + 47re'°'+^pe^~^ / sds 



i.e 



For fc = — 1, 



i.e 



<P + Att f se^^^e^+^pds+ f -L(l _ _^)ds 
Ji 2s^ Co + s' 

< /i + - log(Co + 1) + 47r(Ci + l)2el^l+^p y s'^ds 



p{t) <P+\ log(Co + 1) + 47rel^l+^p(Ci + 1)^ logt. 



We have proven : 
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Theorem 3.2 Consider spatially homogeneous solutions of the plane and hy- 
perbolic symmetric Einstein- Vlasov-scalar field system written in areal coordi- 
nates and initial data given fort — 1. Then these solutions exist on the whole in- 
terval [1,00 [ and the corresponding spacetimes are timelike and null geodesically 
complete in the expanding direction. Estimates f!^.b]} - {^7^ - l^^^.8\} - f^^.y\} - f^^.lU\} - 
l^t^.llp are satisfied. 

Remark 3.3 In the homogeneous plane symmetric case with only a scalar field, 
solving the wave equation, we obtain = ji/'lt From the field equations, 

k{t) = (27r?A2„ i)t-i and \{t) = X+{2n^/P-^)logt ; fi{t) = ti+{2n^lP-\-^)logt 

; Pit) = Pit) = l^j^e-'h-^'^-^' ; jit) = ; qit) = 2p(t). 

We can compute the limiting values of the generalized Kasner exponents : 



lim 



lim 



Kit,r) 
Klit,r) 



*^°o tA + 2 



= lim 

t^oo A (t, r) t^oc 



Klit,r) 
Kit,r) 



= lira 

t — >oo 



3 
1 



tX + 2 4Tnjj^ + 3' 



4 Strong cosmic censorship 

Consider the Einstein-scalar field system with plane symmetry and initial data 
given at t = 1. It follows fi'om 8 that the Kretschmann scalar blows up 
uniformly as the singularity is approached. This means that the inhomogeneous 
spacetime cannot be extended to the past. It cannot be extended to the future 
for any initial data since we have proved geodesic completeness in theorem 13.11 
Therefore the maximal globally hyperbolic development of any initial data is 
inextendible. This means in conclusion that the strong cosmic censorship holds 
for this model (see 4 ). 

The question whether the spacetime is geodesically complete in the non- 
homogeneous hyperbolic case remains open. 
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